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Abstract. —Recent efforts to sustain Burrowing Owl {Athene cunicularia) populations in Manitoba have 
been unsuccessful, and the species is now effectively extirpated from the province. Although specific 
causes of the decline remain unknown, loss, fragmentation, and degradation of suitable habitat have 
likely been major contributors to this decline. We developed a habitat suitability index model to deter¬ 
mine suitability of Burrowing Owl nesting habitat in southwestern Manitoba and southeastern Saskatch¬ 
ewan. Model parameters were obtained using a modified Delphi technique to solicit expert opinions. 

An interactive, adaptive learning approach was used in model development, iteratively refining the 
model until acceptable levels of accuracy and robustness were achieved. Application of the model to 
historical Burrowing Owl breeding sites in Manitoba indicated that habitat suitability is often reduced 
by the presence of tall vegetation at former nest burrows. A management approach involving moderate 
grazing to maintain low vegetation height at all nest burrow sites is recommended. 

Key Words: Burrowing Owl, Athene cunicularia; habitat', habitat suitability index', modeling', grazing, Man¬ 
itoba. 


Desarrollo de un modelo de indice de aptitud del Habitat para los Buhos Cavadores en las praderas 
orientales canadienses 

RESUMEN.— Recientes esfuerzos por sostener las poblaciones del Buho Cavador {Athene cunicularia ) en 
Manitoba no han tenido exito, y la especie esta ahora efectivamente extirpada de la provincia. Aunque 
las causas especificas del declive permanecen sin conocerse, la perdida, fragmentacion y degradacion 
de la aptitud del habitat han probablemente sido los mayores contribuidores a este declinacion. Noso- 
tros desarrollamos un indice de aptitud del habitat para determinar la idoneidad del habitat de ani- 
dacion del Buho Cavador en el sudoeste de Manitoba y el sudeste de Saskatchewan. Los parametros 
del modelo fueron obtenidos usando una tecnica Delphi modificada para solicitar opiniones expertas. 
Un acercamiento interactive, de aprendizaje adaptativo fue usado en el desarrollo del modelo, refinando 
iterativamente el modelo hasta lograr niveles aceptables de exactitud y robustez. La aplicacion del 
modelo a sitios historicos de reproduccion de Buhos Cavadores en Manitoba indico que la aptitud del 
habitat a menudo se reduce por la presencia de vegetacion alta en las antiguas cuevas nido. Se recom- 
ienda un enfoque de manejo que involucre un pastoreo moderado para mantener la altura de la 
vegetacion baja en los antiguos sitios de los nidos cueva. 

[Traduccion de Victor Vanegas y Cesar Marquez] 


Declines of Burrowing Owl {Athene cunicularia) 
populations in Canada (Haug and Oliphant 1990, 
James et al. 1997) have resulted in the species be¬ 
ing designated as endangered (Wellicome and 
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Haug 1995). Population declines have been ac¬ 
companied by a contraction of the Canadian Bur¬ 
rowing Owl range, with the most pronounced 
range reduction occurring in the eastern Canadian 
prairies (Hjertaas 1997). Recently in Manitoba, the 
species has become effectively extirpated (De Smet 
1997, Rothfels et al. 1999). Although many factors 
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could be responsible for these changes, reductions 
in the quality and availability of grassland habitat 
are believed to be the major factors that have im¬ 
pacted Burrowing Owl populations (Zarn 1974, 
Haug and Oliphant 1990, Millsap and Bear 1997). 

Because wildlife-habitat modeling often facili¬ 
tates an improved understanding of the impacts of 
habitat alterations on wildlife populations (Morri¬ 
son et al. 1998), we developed a habitat suitability 
index (HSI) model for Burrowing Owl populations 
in their former range of southwestern Manitoba. 
Habitat suitability index modeling was originally 
developed by the U.S. Fish and Wildlife Service 
(1981) as part of their Habitat Evaluation Proce¬ 
dures. HSI models evaluate habitat in relation to 
environmental factors that are deemed most im¬ 
portant in influencing the presence, distribution, 
and abundance of a given species (Morrison et al. 
1998). Such models can provide a repeatable as¬ 
sessment procedure for identifying changes in hab¬ 
itat suitability over time (U.S. Fish and Wildlife 
Service 1981, Schamberger and O’Neil 1986, Mor¬ 
rison et al. 1998). Our objective was to develop an 
HSI model that could be used to determine wheth¬ 
er grasslands of southwestern Manitoba provide 
the habitat conditions required to sustain Burrow¬ 
ing Owl populations. 

Methods 

Modified Delphi Data Collection. Little information 
has been published on the site-specific, qualitative habitat 
requirements of Burrowing Owls in the eastern Canadian 
prairies. To obtain this information, we used a modifi¬ 
cation of the Delphi technique (Dalkey 1969, Crance 
1987) to solicit opinions from regional Burrowing Owl 
researchers (Crance 1987). The Delphi technique origi¬ 
nally used anonymous questionnaires to obtain informa¬ 
tion from experts and facilitate consensus building, but 
was subsequendy modified to include group discussions 
among experts (Crance 1987). For the purpose of this 
study, the Delphi technique incorporated both a group 
discussion and a questionnaire component. 

Five researchers, who had conducted studies on Bur¬ 
rowing Owls in southwestern Manitoba and southeastern 
Saskatchewan, participated in a workshop on 13 June 
1997 in Regina, Saskatchewan. The workshop began with 
an overview of HSI models and the principles driving 
their development. Participants were then asked to iden¬ 
tify nesting and foraging habitat requirements of Burrow¬ 
ing Owls in southwestern Manitoba and southeastern Sas¬ 
katchewan based on their research experiences. As each 
habitat requirement was identified, participants collec¬ 
tively developed a suitability-inclex (SI) curve displaying 
the relationship between the habitat variable and the in¬ 
dex of suitability. Participants considered each of the 
identified habitat components, and used a secret ballot 
voting method to select 10 components they believed 
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Figure 1. Flow diagram summarizing the development 
of the Burrowing Owl Habitat Suitability Index (HSI) 
model. Curved boxes represent the Delphi expert work¬ 
shop stage, rectangles the numerical model development 
stage. 


were most important in determining Burrowing Owl hab¬ 
itat suitability. Votes were displayed and workshop partic¬ 
ipants then collectively modified these results as required 
to achieve group consensus. 

Participants were then asked to consider interrelation¬ 
ships among the 10 habitat components as they related 
to Burrowing Owl habitat suitability and reproductive 
success. Individually, participants ranked each compo¬ 
nent by considering its importance for suitability of Bur¬ 
rowing Owl habitat. In addition, participants were asked 
to assign relative weights (scaled from 0.0-1.0) to habitat 
components, resulting in a mathematical description of 
the interrelationships functioning within the habitat. 
Equal ranking and weightings among two or more vari¬ 
ables were permitted. Results of this secret ballot were 
subsequently displayed. Participants then reviewed and 
critiqued the component rankings and weights until con¬ 
sensus was reached. 

Results of the workshop were later summarized and 
mailed to the participants for comments (Fig. 1). Partic¬ 
ipants were asked to approve or modify the results and 
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provide clarification where necessary. Comments and 
suggestions from this questionnaire were incorporated 
into the workshop results and again forwarded to partic¬ 
ipants for their final approval. 

Field Data Collection. To examine present-day habitat 
suitability of historical Burrowing Owl nest sites in south¬ 
western Manitoba, fieldwork was undertaken in the vicin¬ 
ity of Mehta (49°10'N, 101°00'W) within the mixed-grass 
prairie ecoregion of Manitoba (Scott 1996). In early-June 
1998, proximate habitat was examined at 13 historical 
nest sites used by Burrowing Owls between 1987-97. Re¬ 
cords of fledging success between 1987-97 (De Smet 
1997, K. De Smet unpubl. data) were used to classify his¬ 
torical nest sites as successful (70% of broods fledged), 
marginally successful (30-50% of broods fledged), or un¬ 
successful (0 broods fledged). Of the 13 sites, six were 
classified as successful, two as marginally successful, and 
five as unsuccessful. The historical nest sites were located 
in present-day cattle pastures subjected to a variety of 
grazing intensities. Sites were interspersed among cereal 
and forage crops, haylands, summer fallow, and other 
grassland habitat types. 

Because most Burrowing Owl activity occurs in habitats 
located <600 m of the nest burrow (Haug and Oliphant 
1990), this study assessed habitat located within a 600-m 
radius of each nest site. Nest and forage vegetation 
heights, burrow availability, topography, perch availability, 
and openness were assessed at each nest site. 

We assumed that Burrowing Owls select forage habitats 
consistent with optimum-foraging theory (Stephens and 
Krebs 1986). Therefore, potential Burrowing Owl forage 
habitat was examined within 600 m of the nest site using 
information provided by the workshop participants. Veg¬ 
etation structure and composition were visually assessed 
by looking outward while walking a 50-m radius circle 
centered on the nest site. Roadside habitat within 600 m 
of the nest site was also assessed visually using informa¬ 
tion provided by workshop participants. Each distinct 
habitat was classified and assigned a suitability index (SI) 
based upon the workshop rankings. The areal extent of 
each habitat class was also estimated. Vegetation sampling 
was then undertaken, beginning with forage habitat hav¬ 
ing the highest SI and proceeding to habitats having low¬ 
er SI values until a total forage habitat area of 9 ha or 
more was achieved for each nest site. Within each habitat 
type, vegetation height was measured at 10-m intervals 
along three randomly-positioned, 100-m transects, and 
mean forage vegetation height was calculated. 

Nesting activity of Burrowing Owls is restricted primar¬ 
ily to habitat located within 50 m of the nest burrow 
(Haug and Oliphant 1990). To identify habitat suitability 
within this area, vegetation height was measured in each 
of the four cardinal directions 1 m from the nest, and 
also at 10-m intervals along a circle with a 10-m radius 
centered on the nest burrow. At all sample locations, 
height of vegetation was recorded to within 1 cm by using 
a meter stick. 

Burrowing Owls on the Canadian prairies rarely exca¬ 
vate their own nesting burrows (Haug et al. 1993),.rely¬ 
ing instead on abandoned badger ( Taxidea taxus) and 
Richardson’s ground squirrel ( Spermophilus richardsonii ) 
burrows (Wellicome and Haug 1995). Participants of the 
modified Delphi process indicated that Burrowing Owls 


typically use holes having an entrance diameter of 8-35 
cm. To determine burrow availability at each site, the 
number of natural and artificial burrows with entrance 
diameters of 8-35 cm were counted within a 10 X 10 m 
random plot located between 10-50 m from the nest bur¬ 
row. In addition, perch availability and habitat openness 
were determined at each nest site by counting the num¬ 
ber of perches and trees, respectively, within a 50-m ra¬ 
dius of the nest. Site topography was assessed visually us¬ 
ing incremental rankings from flat to moderately rolling. 

Model Development. Information obtained from the 
modified Delphi process was used to develop an HSI 
model. To ensure that the identified habitat components 
were appropriate indicators of habitat suitability, the HSI 
compatibility of each variable was assessed using criteria 
developed by Schamberger and O’Neil (1986). In addi¬ 
tion, the SI curves generated by the Delphi process were 
fitted to mathematical functions using regression analysis 
and statistical modeling (Jeffers 1982). 

Model construction. Construction of the HSI model be¬ 
gan by formulating model objectives and assumptions 
Variables considered to be inappropriate to the modeling 
objectives were excluded from consideration. Explorato¬ 
ry data analysis was used to summarize interrelationships 
among the remaining variables, and variables were each 
identified as limiting, cumulative, or compensatory fac¬ 
tors using guidelines established by the U.S. Fish and 
Wildlife Service (1981). A mathematical expression of 
habitat suitability was generated (Fig. 1) using an adap¬ 
tive learning process of combining habitat variables. Both 
arithmetic and geometric means of variables were consid¬ 
ered in formulating habitat suitability. A multiplicative 
application was considered appropriate for limiting vari¬ 
ables, as it ensured that the HSI would equal zero should 
any of the variable SI values equal zero. Cumulative var¬ 
iables were incorporated into the model additively, while 
compensatory variables were incorporated using either 
an arithmetic or geometric mean as deemed appropriate, 
recognizing that geometric means were more sensitive to 
individual low SI values than were arithmetic means. Var¬ 
iable weights identified by workshop participants were ap¬ 
plied to compensatory variables to express their relative 
importance in identifying suitable habitat (U.S. Fish and 
Wildlife Service 1981). 

Model verification. An interactive computer program was 
developed to verify the function, accuracy, and robust¬ 
ness of the HSI model (Fig. 1). Empirical habitat mea¬ 
surements were entered into the program to explore the 
multi-variable behavior of the HSI model, and to ensure 
that the computed HSI values reflected expert opinion. 
If the model behavior was deemed suboptimal, refine¬ 
ments were made iteratively until acceptable levels of ac¬ 
curacy and robustness were achieved. 

Model validation. The HSI model was validated using 
field data from known Burrowing Owl habitat at Moose 
Jaw, Saskatchewan, and from historical Burrowing Owl 
habitat in southwestern Manitoba (Fig. 1). The Saskatch¬ 
ewan data were used to confirm that the model produced 
high HSI values for currently occupied habitat, while the 
Manitoba data were used to assess the suitability of his¬ 
torical Burrowing Owl nest sites for future populations 
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Table 1. Primary habitat variables, variable priorities, 
and variable importance weights (range of possible 
weights: 0.0-1.0), as determined by habitat-modeling 
workshop participants. 


Priority 

Primary Habitat Variable 

Importance 

Weight 

1 

Burrow availability 

1.0 

2 

Forage availability 

1.0 

3 

Vegetation at nest site 

0.8 

4 

Openness 

0.8 

5 

Habitat fragmentation 

0.7 

6 

Forage habitat quality 

0.9 

7 

Inter-nest distance 

0.5 

8 

Areal extent of nest pasture 

0.2 

9 

Topography of nest area 

0.2 

10 

Perch availability 

0.2 


Results 

Model Construction. Workshop participants in¬ 
dividually identified 19 habitat components 
thought to affect the suitability of Burrowing Owl 
breeding habitat in the eastern Canadian prairies. 
Participants then reduced this to 10 components 
thought to be the most important in determining 
Burrowing Owl habitat suitability (Table 1). 

Nine of the 10 habitat components identified by 
workshop participants were specific to habitat 
proximate (<50 m) to nest burrows, while the 
tenth addressed landscape-level habitat fragmen¬ 
tation. Given that the objective of this study was to 
develop an HSI model for small-scale, proximate 
habitat, the fragmentation measure was excluded. 
The remaining nine habitat components were re¬ 
defined as variables for the purpose of model con¬ 
struction. 

The contribution of the nine remaining vari¬ 
ables to the HSI model was determined by consid¬ 
ering the variable interrelationships identified by 
workshop participants. Number of burrows was 
deemed to be a limiting factor in identifying suit¬ 
able nesting habitat because Burrowing Owls rarely 
excavate their own burrow. Because species surviv¬ 
al is dependent upon the availability of suitable 
prey habitat, the quality and availability of foraging 
habitat were also deemed limiting factors. Inter¬ 
nest distance, openness and vegetation height at 
the nest burrow were identified as compensatory 
factors hr cau se high suitability levels for these var¬ 
iables were expected to offset low suitability levels 
of other variables. Topography at the nest site, ar¬ 


eal extent of nest pasture, and perch availability 
were considered to have minimal influence on site 
suitability individually, but collectively, these vari¬ 
ables were determined to be important. Topogra¬ 
phy at the nest site, areal extent of nest pasture, 
and perch availability were therefore considered to 
be cumulative factors. 

An initial model was produced that included all 
nine variables in a weighted geometric mean, using 
the variable weights suggested by the workshop 
participants. The resulting model proved cumber¬ 
some and was insensitive to changes in variable val¬ 
ues. To improve model function, the three vari¬ 
ables weighted lowest by workshop participants 
were removed. Two additional variables were in¬ 
corporated indirectly into the model as stipulations 
guiding model application. Specifically, the HSI 
was automatically set to zero under either of the 
following conditions: 

(1) Openness: tree or tall shrub encroachment 
within 50 m of the nest site. 

(2) Forage availability: no forage habitat located 
within 600 m of the nest site. 

Subsequent model construction and verification 
focused on the four remaining variables. 

HSI models are intended to be general indica¬ 
tors of habitat suitability that are easily and re- 
peatably applied under field conditions. Minimiz¬ 
ing the number of variables in the HSI model 
served two purposes: the model became more eas¬ 
ily applied, and the likelihood of model overfitting 
was reduced (Jeffers 1982). 

Model Verification. Model verification was un¬ 
dertaken by creating an interactive computer pro¬ 
gram to determine individual variable SI values 
and a composite HSI value for a specific set of hab¬ 
itat parameters. Approximately 500 sets of habitat 
parameters were generated to explore the utility 
and robustness of the model. An iterative process 
was used to modify model parameters until a suit¬ 
able model was achieved. The final model took the 
form: 

HSI = [ (S,) (S 2 ) (S 3 0 - 8 ) (S 4 0 - 5 ) ] 1/33 

where S 1-4 are SI values for the following habitat 
variables: 

Sj = burrow availability 

So = foraerr vegetation heiofit 

^ - - ' O O . o 

5 3 = nest vegetation height 

5 4 = inter-nest distance 
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Figure 2. Suitability Index (SI) curves for the four hab¬ 
itat variables used in the Habitat Suitability Index (HSI) 
model. Fitted curves (solid lines) and variable values 
from the Delphi workshop participants (circles) are 
shown. 


Variable weights incorporated into the model 
are provided in Table 1, and suitability index values 
for the four variables are presented graphically in 

Flg - 2 : 

A simpler model, excluding inter-nest distance, 
was developed for nest sites where Burrowing Owls 
were not present: 

HSI = [(S 1 )(S 2 )(S 3 0 - 8 )] 1 / 2 - 8 

Model Validation. Burrowing Owls have nested 
successfully in Moose Jaw, Saskatchewan for a num¬ 
ber of years (E. Wiltse pers. comm.). The suitability 
of breeding habitat was assessed at nest sites locat¬ 
ed within two distinct study areas, the track infield 
of the Moose Jaw Exhibition Grounds and on land 
adjacent to the Lynbrook Golf Course. Eight nest 
sites were examined at the Moose Jaw Exhibition 
Grounds study area while an additional eight nest 
sites were assessed on the Lynbrook study area. 
Burrowing Owls were nesting on both study areas 
when the assessment was conducted. The HSI val¬ 
ues for each of the study areas equaled 1.0, indi¬ 


cating that the model successfully recognized hab¬ 
itat that was being used by Burrowing Owls. 

Of the 13 historical Manitoba nest sites exam¬ 
ined, nine artificial nest burrows and four natural 
nest burrows were used. All historical nest sites 
were located on pastures that varied in the fre¬ 
quency and intensity of cattle grazing. Four histor¬ 
ical sites were being grazed during the sampling 
period, while four others had remained ungrazed 
for a number of years. At the remaining five sites, 
grazing had occurred in the recent past, although 
historical grazing regimes were based solely upon 
supposition. Discussions with landowners indicated 
that these pastures were being subjected to rota¬ 
tional grazing practices at the time sampling oc¬ 
curred. Although other forms of habitat manage¬ 
ment can be used to suppress vegetation height 
(e.g., mowing, prescribed burning), only grazed 
sites were available for this study. 

Individual SI values were determined at each of 
the Manitoba nest sites for three variables: burrow 
availability, vegetation height at the nest, and for¬ 
age vegetation height. Inter-nest distance was not 
included because the sites were not occupied by 
Burrowing Owls when sampling occurred. SI values 
for burrow availability were >0.8 for 11 of the 13 
sites, but two sites had low SI values (<0.5) due to 
a shortage of available burrows. Forage vegetation 
height SI values were all >0.8, and in 10 of the 13 
historical nest sites maximum suitability values of 
1.0 were achieved. SI values for vegetation height 
at the nest showed the greatest variation, ranging 
from 0.19 to 1.0. Together, these results indicate 
that forage habitat quality and burrow availability 
were close to maximum at most sites. By contrast, 
SI values for vegetation height at nest sites were 
often low, particularly in ungrazed pastures. 

HSI values for the 13 sites ranged from 0.58 to 
1.0, and were correlated with historical nest site 
success (r 2 — 0.33, P < 0.05; Fig. 3). HSI values of 
unsuccessful historical sites ranged from 0.58 to 
0.79, with highest values occurring in pastures 
grazed when sampling occurred. By contrast, HSI 
values for successful historical sites ranged from 0.7 
to 1.0, with smallest values occurring in sites that 
had not been grazed for some time. These results 
suggest that moderate cattle grazing of nest sites 
may be critical to the maintenance of suitable Bur¬ 
rowing Owl breeding habitat in southwestern Man¬ 
itoba. 
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Figure 3. The relationship between brood-rearing suc¬ 
cess (1987-97) and Habitat Suitability Index (HSI) values 
(June 1999) for 13 historical Burrowing Owl breeding 
sites in southwestern Manitoba. Success codes are 1 = 
unsuccessful, 2 = moderately successful, and 3 = suc¬ 
cessful. The dashed line was fitted using linear regres¬ 
sion. 


Discussion 

The HSI model developed in this study incor¬ 
porated directly four habitat variables considered 
critical in evaluating the suitability of Burrowing 
Owl habitat in southwestern Manitoba. An addi¬ 
tional two variables were incorporated indirectly as 
stipulations guiding model application. During the 
early model development phase up to nine habitat 
variables were included; however, these more com¬ 
plex models proved to be cumbersome and pro¬ 
duced unstable results. Similar results were ob¬ 
tained by O’Neil et al. (1988). Simpler models 
have the advantages of being tractable mathemat¬ 
ically and easily applied under field conditions. 

Application of the HSI model to historical Bur¬ 
rowing Owl nest sites in Manitoba suggested that 
habitat suitability was most strongly compromised 
by the presence of tall vegetation at the nest bur¬ 
row. In southwestern Manitoba, native mixed-grass 
prairie and tame grasses may exceed 10 cm in 
height by earlyjune, particularly in wet years. By 


contrast, vegetation height in the drier mixed-grass 
prairies of southwestern Saskatchewan rarely ex¬ 
ceeds 5—10 cm, even in the absence of grazing. In 
the Canadian prairies, annual precipitation in¬ 
creases from west to east and results in taller grass 
and forb species in southern Manitoba (Scott 

1996) . 

In Manitoba, moderate grazing appears to be 
critical to maintaining an optimal (<6 cm) vege¬ 
tation height at nest burrows. Historical nest sites 
subjected to catde grazing during habitat sampling 
were identified as having greater habitat suitability 
than ungrazed sites, indicating that grazing may 
enhance the suitability of Burrowing Owl habitat. 
In historically successful breeding habitat, cessa¬ 
tion of grazing resulted in degraded habitats con¬ 
sisting of tall and lush grasses. De Smet (1997) 
identified over 700 pastures in southwestern Man¬ 
itoba as potentially suitable habitat for Burrowing 
Owls; however, the quality of these habitats has not 
been formally assessed. The HSI model developed 
in this study can be used to assess these pastures, 
and identify the variables that compromise habitat 
suitability. 

Population declines in Manitoba have been at¬ 
tributed previously to vehicular mortality and re¬ 
cent inclement, wet spring weather (De Smet 

1997) . Based on our analysis, we suggest that hab¬ 
itat deterioration at historical nest sites may have 
contributed to the effective extirpation of Burrow¬ 
ing Owls from Manitoba. Vegetation management 
at historical nest sites, such as regular grazing or 
mowing and either mechanical or fire manage¬ 
ment directed toward the removal of encroaching 
woody plants, should be implemented to improve 
Burrowing Owl habitat by maintaining consistently 
low vegetation height at nest burrow sites. 
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